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1. Introduction 

It is well knowiiEI that instantons do spoil the operator product expansion by 
introducing corrections that are power-suppressed by 9 or more inverse powers of 
the momentum. Such corrections are fully calculable. In view of the large power 
suppression associated with them, in practical applications they will turn out to be 
either very large, or very small. They provide therefore a lower limit to the momentum 
scale at which perturbation theory can be applied. In a previous workB, the effect 
of small-size instantons on the correlator of two currents was computed, and the 
calculation was applied to the hadronic decays of the r lepton. In this case it was 
found that the effect was small, mainly because of a chiral suppression factor. In 
a subsequent workB it was instead shown that in the sum rules that are usually 
employed to determine the mass of the light quarks, instantons effects are indeed 
very large. 

The aim of the present work is to extend and complete the work of ref. 0, with the 
inclusion of corrections to the coefficient function of quark condensates in the operator 
product expansion of two currents. Such corrections are less chiral suppressed than 
the one one computed in ref. |0] (a discussion of this issue can be found in ref. [??]). 
Results in this direction have already been obtained in ref. 0, where the instanton 
correction of the coefficient function of a six-quark operator was computed in the 
factorization approximation. 

Our paper is organized as follows: in section 2 we discuss the general framework for 
computing instanton corrections to the coefficients of the operator product expansion 
of a current-current correlator. In section 3 we describe in detail the computation 
of the instanton corrections to the coefficient function of the chiral condensate. In 
section 4 we extend our calculation to four- and six-quark operators. In section 5 we 
collect the final results. In section 6 we discuss some phenomenological applications, 
and in section 7 we give our conclusions. 

2. Instanton Corrections to the 
Current- Current Two-Point Functions 

We begin by establishing a framework for the computation of instanton corrections 
to the coefficient functions of the OPE. Let us consider the time ordered product of 
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two currents. The operator product expansion reads (all indices were dropped for 
simplicity) 

J{x)J{y) ^ 5:C«(A)0,(r) (2.1) 

i 

where 

r = -^—, ^ = x -y, 

Oi = 1, O2,... = operators of higher dimension. (2-2) 

In order to compute the coefficients of the OPE in perturbation theory, one usually 
computes the expectation value of both sides of eq. ( |2.1| ) among suitably chosen 
perturbative states. Alternatively, one can compute in perturbation theory both 
sides of the equation 

(0| J{x)J{y) P |0) ^5:C«(A) (0| 0,(r) P |0) (2.3) 

i 

where P indicates a generic product of local fields at fixed positions away from the 
origin, and then solve for the Cq . 

In order to determine the effect of instantons, we may now compute both sides of 
eq. (|2.3|) by including instanton effects. We must have 

(J|J(x) J(i/)P|J) ^EC^^^HA) {I\0.{r)P\I) (2.4) 

i 

where now |/) denotes the vacuum in the presence of instantons, and C*^*^ = Cq +(5C*-*^ 
represent the modified coefficient function in the presence of instantons. In the dilute 
gas approximation we have 

(/|J(x)J(y)P|/) = (0|J(x)J(y)P|0) + 
diD{i) [(e|J(x)J(y)P|O-(0|J(a;)J(y)P|0)] 

(2.5) 

where |^) represents the one-instanton configuration. The variable ^ represents here 
all degrees of freedom (position, size, orientation, colour, signature) of the instanton 

I di= I dp<fzdRY, (2.6) 
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and D{C,) is the instanton density. For the operators we also have 

{I\Oi{r)P\I) = {0\Oi{r)P\0) + 
'd^DiO mO.{r)P\0 - {Omr)P\0)] . 



(2.7) 



From eqs. 
we get 



pTSP and ( p.7| ), treating the instanton correction as a small correction, 



^5c«(A)(o|a,w^|o) 



d^ D{0 



(e|J(x)J(y)P|0-ECo (A)(e|0.(r)P|O 



(2i 



The coefficient functions are then obtained by appropriately choosing P. If we choose 
P = 1 (the identity operator), then only the identity operator O^o will survive on the 
left-hand side, and we will get 



5C^'\A) = I d^ DiO 



(e|J(x)j(y)io-E^o (A)(e|0dO 



(2.9) 



The second term in the square bracket of eq. (|2.9| ) subtracts the effect of large-size 
instantons from the ffist term. In fact, the expectation value of the product of two 
currents in a slowly varying classical field obeys the operator product expansion 



(2.10) 



(e|J(x)j(y)io^E^o(A)(e|0dO 

when the x — > y limit is taken at fixed instanton size. Therefore, in eq. ( |2.9| ) the 
subtraction effectively cuts off the effects of large-size instantons. By looking back 



to the derivation of eq. (2^), it is easy to see that large instantons affect the matrix 
elements of the OPE, but not the coefficient functions. 



Equation ( p.9| ) cannot be valid for operators of arbitrarily high dimension. In fact, 
after integrating over the instanton position, the expectation value of the operator 
Oi in the instanton background is given by an integral over the instanton size, which 
behaves like /c/pp®+"^/'^p~*, where di is the canonical dimension of the operator. 
Therefore, ii di > 7 + ^//3 the integral is ultraviolet divergent. One should then 
worry about subtractions. If we limit ourselves to operators with dimension smaller 
than 7 + nf/3, no UV divergences are present, and the instanton correction to the 
coefficient function is well defined, because all infrared divergences in eq. ( p.9|) are 
removed by the subtraction term. In the particular case in which nj is an odd multiple 
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of 3, not all the infrared divergences are removed in this way. In fact, operators 
of dimension 7 + ?t-//3 may exist, and their matrix elements are both infrared and 
ultraviolet divergent. One then has an ambiguity in deciding the scale of separation 
of short and long distance effects in the OPE. For our purposes, this problem will turn 
out to be irrelevant. In fact, we will see that the ambiguous term is a polynomial in A, 
and therefore its Fourier transform is concentrated at p = 0, and it never contributes 
to physical quantities. 

3. Instanton corrections to the coefficient function 
of the chiral condensate 

In order to compute the correction to the coefficient function of the chiral conden- 
sate C(//) = 4'f{t)ipf{t), we will choose P = ■ipf(t)ipf{t), where the point t is chosen 
to be far away from the origin with respect to x, y. We get 



y: ci^'\^mo.ir)^fit)m\o 



(3.1) 



It is easy to convince oneself that the only operator that can appear on the left- 
hand side is C?(//)- In fact, since we are working at the tree level, in order to get a 
contribution the operator must be a fermion bilinear. Bilinears of the form ipf'-jaipf, 
'i/jfacrs'ipf, etc., give zero after one takes the fermion trace, while bilinears of the form 
ipfd^ipf vanish at a faster rate for large t. 

We will now focus on the computation of the first term inside the square bracket 



of eq. ( |3.1| ) . We define 

i,. = {^\J,ix)JAy)Mt)^ijfit)\0 (3.2) 

where our currents will be in general axial or vector currents, possibly fiavour non- 
diagonal. The Euclidean fermionic propagator in the instanton background satisfies 
the equation 

(z7 -D^ -m) S^{x, y) = 5\x - y) . (3.3) 

In the m —^ limit it has the expansion 

S^(:,^y) = ^tMMyl + s^{x,y) + m f d'z S^{x,z) S^{z,y) + 0{m') ; (3.4) 

lit J 
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5*0 satisfies the equation 

where D^ is the covariant derivative in the instanton background 

In operator notation, 5*0 has the exphcit expression^ 

/ _^2 2 -L)2 ' 2 ' 

where the meaning of the operator notation is specified as follows 

d^xd'^yS^{x,y) /i(x) f2{y) = J d^x d'^y fi{x) S^f2{y) . 

The null-mode projector 

k 



is given in operator notation by the formula 



7-D 



-D'- 



■7-L) 



1^75 



Another useful identity is the following 



m 



n 1 „ 1 1±75 1 lT75 



-D^ ' -D^ 2 -D^ 2 ■ 

Equations ( p.7|) , ( p.lO|) and ( |3.11| ) are easily verified, using the identity 



, O 1 ± 75 r, 1 ± 75 



D' 



2 2 

which is valid for any self-dual (anti-self-dual) field configuration (see ref. 

In operator notation the fermion propagator is then 



S^ = + ^0^ + m (S^Y + 0{ 



m 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



!)• 



(3.13) 
Let us assume now for concreteness that our currents are ud currents, and that the 



scalar bilinear is a dd current. We have 



'■/lU 



-Tr 



T,{x\S^\y)T,{y\Sm{t\S^\x) +Tr T,{x\SM^^{y\S^\x) 



xTr 



{t\sm 

(3.14) 
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(the minus sign in the first term is due to the fermion loop). Observe that we do not 
need to include tadpole diagrams involving the currents, since we are considering the 
flavour off-diagonal case. We want to isolate the least chiral suppressed contribution. 
We immediately find that the relevant contributions are given by 



'■/ii/ 



nd ^ mi ^ m„. 



rrid 



r' 



(3.15) 



with 






fib) 

Tie) 



Tr 
Tr 
-Tr 



V,{x\Si\y)VMS^\t){t\P\x) 
V,{x\S^\y)V,{y\P\t){t\S^\x) 

TM{s^f\y)^MP\t){t\p\x) 



Tr 



TMP\y)TMSm{t\Si\x) 



The remaining term is the tadpole 






Tr 



vMst\y)TAy\sf\x) 



X Tr 



{t\si\t) 



(3.16) 



(3.17) 



which is proportional to the vacuum term. All other terms are either more suppressed 
by powers of the fermion masses, or they vanish because of a wrong chiral structure. 

The computation was carried out in the regular gauge. This is legitimate, since 
I^j_y is a gauge-invariant quantity. The instanton field in the regular gauge has the 



expression 



A^{x)= AUx)t' 



1 V, 



fii/a 



Xv (^n 



g x^ + p^ 



(3.18) 



where cxa are the three Pauh matrices in the SU{2) subgroup of S't/(3) in which the 



instanton lives. Equation ( p.l8| ) is for an instanton centred at the origin (we can go 
back to the general case by the replacement x^x — z,y-^y — z). The 77 symbols 
are defined in ref. [0. The scalar propagator has the simple expression 

1 , , 



{x\ 



-D^ 



P^ + x-y + ir]^^^ Xf, y^ a^ 
An"^ (x — yy (x^ + p^)2 (y^ _)_ p2^2 



We introduce the r symbols 



T„ 



10-1,0-2, org, 2j 



(3.19) 

(3.20) 



The rj symbols satisfy the following equations 



tItu 



^P^^l 



^l^y + ^ Tlt^ua ^a 



^i^v + i Vt^ua ^a • 



(3.21) 
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In terms of the r symbols we have 
, , 1 , > 



[x\ 



p2 + (rt ■ x) (r • y) 



-D2 



47r^ (x — yY (x^ + p^) 2 (y^ _)_ p2^ 2 



i'4"(^) 



% X ^ T^ T ' X 

g x^ + p^ 



(3.22) 



(3.23) 



(y4^(x) e A^{x,y) are obtained with the substitution r ^^ t^). It is also convenient 
to find a simpler expression for the null-mode projector P. The zero modes in the 
regular gauge are given by 



^ [X^ + p2] 2 



Z74 72 



1-75 



(3.24) 



where a is a spinor and i is a colour index. The index i spans the two-dimensional 
space of left-handed spinor, and it corresponds in colour space to the SU{2) subspace 
in which the instanton lives. To be more specific, let us use the following representa- 
tion for the gamma matrices 



so that 



We have 



Therefore 



7 



r 



-i 
-i 



7 







I'Y 



r^t 



-r'^V'^ 



7 



r^ 
-f 

-1 
1 



-T-M^^^t 



r^V*^ 



1 -75 



YY 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



where the left-hand side is naturally extended to four dimensions, with all the entries 
equal to zero, except for the upper-left two-dimensional block. 

Formula ( p. 24] ) was obtained via a gauge transformation of the standard singular 



gauge expression of ref . [T| 



p 



1 

^ [x2 + p 



21 f 



Z7 ■ X 1 

~\ — r 72 — 

X 



75 



With the above normalization we get0 



d X Tr 



ipo{x)^l{x) 



(3.29) 



(3.30) 



^The present work differs from ref. [|l4| in the normalization of the zero mode. 
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and the null-mode projector is given by 

{x\P\y) = ilJo{x)i/jl{y). 



(3.31) 



Consider now a generic diagram containing a fermionic loop. If a null-mode pro- 
jector is present, colour and spinor indices will mix, and Dirac and colour traces will 
generally have the form 



= [7ai • • • lajf3a [Mx)U [4^(3, ■ ■ ■ rl^^jpjjj [My)]](3 , 

Using now eq. (|3.28|) , together with the transposition rule 



T 

7m 



- 72747^7472 , 
we can transform our espression into a standard gamma matrix trace: 

■1-75 



(3.32) 



(3.33) 



-1 - Tr 



P^ 



(7ai • • • 7aJ Mx) 72 74 hfSm • • • 7/3i) 74 72 -iPoiy) 
1-75 



7r2(x^ + p^) 2 (y2 _)_ p2-) 2 



:-l)- Tr 



■ (7ai •• •7«„) (7/3™ •• -7/31) 



(3.34) 



Notice that the factors 7274 in this formula cancel against the analogous factors 
appearing in the zero-modes, thus yielding a fully Lorentz invariant result. The spin 
and colour traces one encounters when computing the quantities ij^j and ijfj are 
precisely of the above form. In the case when we have two null-mode projectors in 
the trace, as in the ijj'j term, we get instead 



(7ai • • •7an)/3a i^Jp^ . ..rl^_jpjji [Mx)U [Mt)]*Sk[Mt)]sk[My)]}j 



2p' 



7r4(t2+p2)3(a;2+p2)f(^2 + p2)f 



Tr 



1-75 



(7ai---7a„)(7/3„---7/3i) 



(3.35) 



With all the machinery developed so far we can perform the calculation using standard 
algebraic tools. We need to perform the spin and colour trace first. We choose 
for convenience a frame in which x = —y = A/2. In the approximation |t| >> 
|x|, \y\, \z\,p, neglecting terms that have opposite sign for the instanton and for the 
anti-instanton, we get 






p^ 



2 7r6t6A4 



1 1 

+ 



dA.d 



dPid. 



A^6 



flU 



A„A, 



2A^ 
dldl 



■A.A, 



(3.36) 
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Tic) 



^2 



P 



-V^^e^eV 




(3.37) 
(3.38) 



d-i 



2 
and ip is defined to be 1 for vector- vector and 



where we have defined 

+ p' (3.39) 

1 for axial- axial currents. We do not 
need to take into account the region of integration where 1^1 is of the order of |t|, or 
where 1^ — t| is of the order of p. According to our formula ( p.l|) , such corrections 
will be subtracted away, since in fact they can be viewed as instanton corrections 
to the expectation value of the operator, and not to the coefficient functions. The 
integration over the instanton position can be performed as follows. We define the 
integrals 

1 „ ./•,.! 



d'^z- 



H 



d^'z- 



d+d- ' J d+d"^ 

The following identities are easily proved 

OF „ Id 



G 



d'^z 



1 



dldl 



H 



G = - 



A' 



dF 



(3.40) 



(3.41) 



Qp2' p2Q^2 y Q^2j ' 

Observe that F is logarithmically infrared divergent, while H and G are finite. Since 
F is dimensionless, it will always be possible to express it as a function of the ratio 
A/p plus a term proportional to log -^, where L is the infrared cutoff. On the other 
hand we can easily verify that both H and G remain invariant if we add to F a term 
proportional to log A. Therefore we can always choose F as a function of the ratio 
A/p alone. If we restrict ourselves to this choice we also have 

dF A2 dF 



H 



d p^ 



p2 9A2' 



(3.42) 



We have therefore 



Tib) 



Tic) 
fj,u 




(3.43) 



'"AtiH^^"" 
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We only need to compute F. We get 



F = 7rM-log4^ + iH 



e + 1 



(3.44) 



with ^ = wl + AlP^ 1 1^. The Melhn transform of F is 

-2 ^ M+l\ . /M+1 



"M 



p^'^Frfp = -A^'^+^ — COS vr 



A/+1 



vr 



r^ 



r(-M-2)(M + l). 

(3.45) 



Using the formula 

/'d^Ae'^PA^=p-"-^47r 



. 7r(iv + 2) r(f)r 

sm — ^-^ — 



iv+2^r(Ar + 4) 



r(^) 



(3.46) 



it is now easy to perform the Mellin transform and the Fourier transform of our result 
by a simple algebraic procedure. Defining 



we get 



i(;M(M) = ld'Ae^^-^ Jdpp'^li^/ 



lSiM) = -^p-''-'p,P. 



b,c) 



iSiM) 



MM) _M-7 



\t^7T^ 



P 



5^vP^ 



~mM)=^,^p-''-'S,.p' 



where 



A{M) = 71^ 



r(|)r3(^)(M + 3)2(M + 5) 



A/+6 



(3.47) 

(3.48) 
(3.49) 
(3.50) 

(3.51) 



The calculation of the disconnected term (eq. (|3.17|) ) is quite simple. We have 



Tr 



{t\sm 



.^Tr 



{t\P^\t) 



2p2 



md TT^ t^ 



(3.52) 



The remaining factor \r ^{x\S^\y)V u{y\S^\x)] can be easily obtained by suitably 
adapting formula (2.28) of ref. 0. We get 



C(M) 



A{M) p 



-M-7 



AnidirH^iM + 6) 



M + 5) ipC 5f,yp^ - p^p^ 



+ {ipC-l)5^,p'], 

(3.53) 
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where C = {iriu/md + md/mu)/2. The relative minus sign with respect to ref. P| 
is due to the fact that in our case the minus sign of the fermion loop has cancelled 
against the minus sign of the fermion loop of the tadpole, and the extra factor of 1/2 
accounts for the fact that our expression is not yet summed over the instanton and 
anti-instanton. We now need to include the instanton density in our result. We follow 
closely the notation of eq. (3.10) of ref. 0, defining 



D{p) = H 



log 



p2A2 



P 



6+-f 



(3.54) 



where H and c are given in ref. Pj. In leading logarithmic accuracy (which is the 
accuracy of our calculation) the logarithmic factor in the integration can be taken 
out of the integral sign, by simply replacing p with 1/p. We get 



di D{i) /, 



^v 



A(y)p-^^ 



H\\og 



p- 



it^n^ 



rrid 



livM-'. 



Pf^Pu + ip 



1 



-^ + — ]p s^ 

mi m„. 



/ii/ 



i^uP^ - Pt^Pv) + -^{hC -I) 6^^p' 



. (3.55) 



Observe that in the case m„ = rrid the result is transverse, which is an important 
check of our calculation. Another important observation has to do with the relative 
sign of the tadpole term with respect to the rest of the expression. Observe in fact 
that there is a term proportional to mu/m'^ in the expression, which cancels after 
inclusion of the tadpole term. Since H contains the factor mumdms, such a term 
would otherwise give rise to corrections proportional to rris mu/md, so that the chiral 
limit rud, rriu — * (with nis fixed) would be undefined, with disastrous consequences 
for the usual interpretation of chiral symmetry in QCD. 

Observe that in intermediate steps of the calculation there are infrared divergences 
that do not appear in the final answer. In fact, the Mellin transform of F, eq. (|3.45| ) 
is divergent for odd integer M. In particular, for three light flavours M turns out to 
be an odd integer. Our final answer is however finite, because the Fourier transform, 
eq. (|3.46| ), vanishes when N is an even integer (in fact, the Fourier transform of an 
even power of A is a derivative of the four- dimensional delta-function of p, and it 
vanishes when p is away from zero). One may therefore doubt that our result may 
depend upon the regularization method that we implicitely used, which is to continue 
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our result for non-integer M. We should however remember that the IR divergent 
terms should be subtracted from I^^, according to eq. (|2.8|). As already discussed in 
section 2, the subtractions are defined up to finite terms which are polynomials in 
A. Once the subtractions are performed, the Mellin transform will turn out to be 
infrared finite, and it can therefore be regulated in any way we like. We can compute 
it for complex M, perform the Fourier transform, and then take the appropriate limit 
for M integer. If we use this procedure the subtraction terms do not survive, because 
they are polynomials in A, and therefore they have zero Fourier transform (for a more 
detailed discussion see ref. [@]). 

Following eq. (|3.1| ), in order to give our final expression for the correction to SC^'^'^'^ 
we still need to divide by 



it ■ 7 —it ■ 7 



27r2 14 27r2 f^ 



nH^ 



(3.56) 



(0|O(,-,)(0)^,(t)^,(t)|0) = -3Tr 
(the —3 comes from the fermion loop and the colour sum). Our final result is then 

5^^) = ^^ (\ok^Yp~''— {6,.p'-P,P.) . (3.57) 



flU 



rrid 



The corresponding expression for SC^^'^^ is obtained with the obvious replacement 
nid — ^ rriu. The coefficient SC^^'^^ instead receives contributions only from the tadpole 
term. We have 

-2 ', 

Pf,Pu) + (ip C 



5Cf,ff^ 



flU 



HA{7) 
' 787r2 



-F, 



rus 



12UpC6^,p' 



l)5,.up' 



(3.5J 



The case of flavour diagonal currents can be treated similarly. For a dd vector current 
we get 

P'^^ — fV/ -PtiPu + 13 {ipP^S^y - p^Pu)] (3.59) 



_ HA{1) r 



y{dd) 



rrid 



liUjdd 



HA{7) 

■ 787r2 



log 



p^ 

A2 



p 
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^,iuP'^ -Pt^Pu + 13 (ipP^S^y 
12 [ipd^up^ -PtJiPu) + {i 



i)Vp' 



(3.60) 
and a similar one for the uu contributions. Observe that in the case of axial currents 
eqs. ( p.59| ) and ( p. 601 ) are incomplete, and other contributions must be added. This 
is because in the flavour diagonal case we could have extra diagrams, in which there 
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is a tadpole attached to one or both of the currents. It is easy to show that such 
contributions vanish for the vector currents. In fact, given the tadpole expression 

T, = Tr[T,{x\S^\x)] (3.61) 

rotational invariance implies that it must have the form (taking the instanton centre 
at the origin) 

T^ = x^g{x^). (3.62) 

In the case of a vector current we must also have 5^T^ = 0, which implies immediately 
T^ = 0, since we cannot have a radial field with no sources. Therefore there are no 
extra tadpole contributions in the case of vector currents. The divergence of the axial 
current is instead non-zero, and its value is determined by the anomaly. Since we 
have no application in mind for the flavour-diagonal axial-current correlator, we will 
not extend our calculation in order to cover this case. 

4. Instanton corrections for higher-dimension condensates 

The coefficient functions of operators with four and six quark fields also receive 
corrections from the instanton. In particular, the coefficient functions of four quark 
operators can receive corrections of the order of rris, while the coefficient functions 
of six-quark operators can receive corrections with no chiral suppression factors at 
all. Since there are many such operators, and since their expectation values are not 
known, it would be useless to determine the corresponding correction to each coeffi- 
cient function. In ref. [^ the simplifying assumption was made that the expectation 
value of six quark operators factorizes in terms of the J2f '^f'^f quark condensate. 
Under this assumption, we do not need to determine the corrections to each coeffi- 
cient function independently. In general, in order to compute the effect of 2nf quark 
condensates in the factorization approximation it is enough to compute eq. (|2.8| ) with 
P = {J2f V'/^/)"^ iu order to generate on the left-hand side of the equation the com- 
bination of operators appropriate to the factorization hypothesis. In other words, 
we use the operator P as the source of a factorized expectation value for 2 nj quark 
operators. We then observe that 

(0|(nA(0)^/(0)) (E,^/W^/W)"^|0) = nj\ (^-^y (4.1) 

while on the physical vacuum, according to the factorization hypothesis 

(vac|n,^/(0)V^/(0)|vac) = (gg)"/, (4.2) 
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where we have assumed a flavour-symmetric vacuum expectation value 

(vac|^/(0)^/(0)|vac) = (gg) (4.3) 

independent of /. In order to get the correction to the correlator coming from n/ 
quark condensates in the factorization approximation we should therefore multiply 
our result by 

In order to simplify the calculation, one should keep in mind that terms proportional 
to the inverse of a quark mass cannot appear with a power greater than 1, since they 
must cancel against the factor of a quark mass in the instanton density. This fact is 
rather general, and it has to do with the fermionic nature of the zero modes. The 
contributing terms are given by the tadpole-type term 

Kl^J = -TT[T,{x\S^\y)^T,{y\S^\x),] x Tr[(t|P|t)J Tr[(t|P|t),] Tr[(t|P|t),] 

(4.5) 
(where the 6 is a combinatoric factor) and the following connected contributions 

Kf} = —TT[{t\P\x),r,{x\S^\y)^r,{y\S^\t),] x Tr[(t|P|t)J Tr[(t|P|t). 

+ —Ti[{t\S^\x)dTMS^\y)u^Ay\P\t)d] x — ^Tr[(t|P|t)„]Tr[(t|P|t), 
nid rriums 

+ 1. TT[{t\P\y)^r,{y\S^\x)ar^{x\S^\t)^] x — ^ Tr[(t|P|t)d Tr[(t|P|t), 
m„ rridms 

+ — TT[{t\S^\y)^r,{y\S^\x),r^{x\P\t)u] x — ^ Tr[(t|P|t)d Tr[(t|P|t), 
rriu mams 

K\^ = —TT[{t\S^\x),T,{x\P\y)^T,{y\S^\t),] x -^ TrKt|P|t),] Tr[(t|P|t), 

+ —TT[{t\S^\y)uT,{y\P\x)dT^.{x\S^\t)u] x — ^ Tr[(t|P|t)J Tr[(t|P|t), 
rrid rridms 

KlfJ = —n{t\S^\x),T,{x\P\yUTAy\S^\t),{t\P\t),]x^TT[{t\P\t),] 

2 „ ,,., ±1 , /„,iDi^\ -p /^ic±i-> ' ' ^ 3 



TT[{t\S^\y)uT,{y\P\x)aT,{x\S^\t)^{t\P\t)^] x — Tr[(t|P|t), 
Ki'J = —TT[{t\S^\x),T,{x\P\tUt\S^\y)uTAy\P\t)d]y<—TT[{t\P\t). 
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—TT[{t\P\x)aT^{x\S^\tUt\P\y)uT,{y\S^\t)d] x A Tr[(t|P|t),] 

niumd nis 

(4.6) 

The suffixes u, d, s are shown as a reminder of the flavour flowing in the line. When a 
projector P is present in the spin-colour trace, the trace can be split into expectation 
values over zero modes, by replacing P -^ ipip'^ . Defining as usual 

k\^-'\M) =Jd'Ae^^-P Jdpp^'Kl^-^^ (4.7) 



we get 






m 



A{M')p-^'-' ( ,^„ ^. , 2\ 6 



Si#Ti){-(^''' + **-*-^^}^^ 



morris Vvr 



A{M').,_,_ ^^ 12 



P PfiPu X 



F 4r°7r° rriumdms \ti 

f>ic) _ ■ MM') -M'-7 X „2 



A(M') _M'-7. 2.. 12 / 2 ^2 






K^H} 



A(M') _^,,_7 12 / 2 



-*p^7?ir?rP PmP-x 



2 



where now M' = M + 4 (because of the four extra powers of p coming from the 
tadpole terms). The term KJ^J was obtained from eq. ( p.53| ) suppressing the term 
proportional to C, while KJI',} and KJ^J are taken from eqs. (|3.48|) and (|3.50|) . Only the 
terms Kjfj and KJ^J have a form different from that of the previously computed terms. 
Their calculation can however be easily performed using the techniques previously 
described. Adding up all the contributions we get 

.i.e (m„mdm,) Trio ti8 M' + 6 

> + ^i^p - 1)] {P,.P. - W') • (4-9) 

Taking now into account the normalization factor eq. ( |4.4| ) and the instanton den- 
sity, including a factor of 2 for the instanton anti-instanton contributions, we get 
the correction to the vacuum polarization coming from six quark operators in the 
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factorization approximation 

'''' \ A^J 51 murridms 

X [1 + i|(2p - 1)] (p^p, - (5^,p2^ . (4.10) 

We observe that with no further effort we can also write down at this point the 
correction to the vacuum polarization coming from four-quark operators involving u 
and d quarks only in the factorization approximation. This is a sensible thing to 
consider, since the strange quark mass is not so small, and one may doubt that a 
factor of l/nis can provide a sensible chiral enhancement. We get 

'^^ X^k^) 135 m^rrid 

x[l + i|(2p-l)](p^p,-(5^,/j_ (4^^^) 



5. Final results 

It is now time to collect all our results. We should remember that the quark 
masses appearing in our formula are running masses, given by 



m(/i) = m ( log — j . (5.1) 

The condensate is also renormalization-group dependent. We have 



(gg) = -A^(log^)^^-^"^ (5.2) 

In terms of the renormalization-group-invariant quantities, we get 

. 24 

x< [ — + — + —]Uf,yP^ -p^py\ - — {ipCdf^^p^ -p^Py\ \ (5.3) 

. 48 



f"" 135 V ^y \P^ P' J 

x[l + i|(2p-l)](p^p,-Vp') (5.4) 
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51 



P 



2 \ ^^+5 



9 r,9' 



2 — / log^ 



'"/ / A^ ft 



p9 p9 



[1 + Ti(^p-1)] {Pl^Pu-^i^uP^ 



(5.5) 



where 



i^n = ^ 



i( ^^ j 2««--f exp[-a(l) + i + (2nf-2)a(i)] 



TT 

45 - 5nf 



33 - 2n 



C 



1 /^ ^\ 



9 V + / 

2 Vm^ niuJ 

1 for vector currents, — 1 for axial currents. 



(5.6) 

—g^u in order to 

go to the Minkowski space, our formula for the current- current correlator, including 
instanton corrections is then 



For Uf = 3, performing the replacements p"^ 



-p^ and 5^ 



flU 



^f.uip'' 



Atx' 



2 M / 2n m„mdms / 5.1701 a 
[Pi^Pu-p 9i.vnog[-p)- 



9 r 



po 



P 



log 



A^ 



[%ipC - To) iP^9t^u - Pt^Pu) + ^,{ipC - l)Pf,Pv 



ThuThdrhs fi^ /6.3718A^ 



po 



p 



log 



—p^ 



-\2 



— + — + —) (Pt^Pu- 9^,vP^) + ^— (ipCg^^p^ -p^Pu) 



m„ mrf m. 



m./i'' / 14.446 A 



pi 



P 



P? p5.370A 



P' 



P 



log 



log 



— p^ 



-p^ 
A^ 



m.. 



a + i(v 



1)) (p. 



Vv- 9i,vP 



■) 



;i + Ti(^p- 1)) [Pf^Pu-g^uP 



(5.7) 



where we have also included for completeness the term computed in ref. [0]. The 
last term of our result agrees with ref. P], except that the expression they use for 
the instanton density, taken from ref. ^, does not agree with ours. The instanton 



density for SU{N) was computed in ref. |10| in the Pauli-Villaars scheme, and the 
conversion to the MS scheme was performed in refs. [Q and |Tl| (for a more detailed 
discussion see also ref. 0). 
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In order to estimate the instanton corrections, we choose the following values for 
the various parameters 



m„ = 8.7 MeV, m^ = 15.4 MeV, m, = 283 MeV, 
fi = 180 MeV, A3 = 400 MeV. 



(5.^ 



We find then 



47r2 



{P^iPu - P^g^^u) log(-p^) - 



0.738 



12 



log 



-p" 
A^ 



X 



{% hC - To) {P^a^iu - Pi^Pu) + ^>{ipC - l)p^,Pu 



'0.886l\ 
. P ) 



14 



log 



-p^ 



A2 



X 



Pt^Pu- 9i,vP ) + 



0.25f2„C(7^,/ 



Pt^Pv 




log 


A2 _ 


log 


VI 

A2 _ 



16 Up 



M 



Pl^Pu-QtiuP 



(l + f(ip- 1)) [P^^Pu-9f,uP 



(5.9) 



This result deserves a few comments. First of all, we see that the six-quark correction 
is of the same order of magnitude as the four-quark correction, they both become 
of order 1 at a scale of 1.3 GeV. The two-quark correction becomes of order 1 at a 
scale of 0.9 GeV, while the vacuum correction is important at a lower scale of 0.74. 
This is consistent with the assumption that QCD is very near the SU{2) chiral limit, 
but only marginally near the SU{3) chiral limit. Therefore corrections suppressed by 
powers of m„ or rria are indeed small, while those suppressed by powers of rris are not 
necessarily small. We also observe that although the four- and six-quark corrections 
are larger than the vacuum and two-quark corrections, they do not always dominate. 
In fact they are purely transverse, and they do not enter in sum rules involving the 
divergence of the current. 

The four- and six-quark corrections have been obtained by assuming factorization, 
while this assumption was not needed for the two-quark correction. We should there- 
fore consider the latter as being on more solid ground than the former. Nevertheless, 
it is fair to assume that the corrections computed in the factorization approximation 
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should at least give an order-of-magnitude estimate of the effect. We also notice that 
the six quark correction has the same power behaviour of the two instanton correction 
(see ref. g). 

In ref. [^ the integration over the instanton density is performed with a method 
that includes also effects subleading by powers of logarithms of p^/A^. In practice 
these effects lower the effective scale at which the logarithms are computed. We 
preferred instead to use only the leading logarithmic expression. In fact, if p ^ A one 
may indeed get a better estimate of the effect by using the method of ref. . However, 
in the case when the effective scale is too close to A, the logarithm approaches zero, 
and it therefore yields a suppression that does not have any physical basis. In other 
words, we should regard our logarithmic factors as expression that do in fact approach 
a logarithm for large values of their argument, but become of order 1 when the 
argument is of order 1. 

In the case oi V — A currents, the terms proportional to ip disappear, and in the 
four-quark and six-quark corrections a further suppression arises, which can be viewed 
as a cancellation between the axial and vector current effects. As noticed in ref. [H 
such cancellation could be an artefact of the factorization approximation. In the 
present work we prefer to take the results in the factorization approximation at their 
face value, in order to avoid making too many uncontrolled assumptions. Similarly, 
we completely neglect the fact that the factorization hypothesis is inconsistent with 
the renormalization group equation. 

6. Phenomenological results 

In this section we will discuss some phenomenological applications of our calcu- 
lation. First of all, we would like to stress the basic difference between the formulae 
for instanton corrections and standard perturbative formulae. Instanton corrections 
are proportional to a power of A. Therefore, at the leading logarithmic level they 
are defined up to a multiplicative constant of order 1. At the next-to-leading or- 
der level, the prefactor is defined only in a leading logarithmic sense. The powers 

2 

of log j2 appearing in the prefactor can therefore be written in various ways, differ- 
ing by subleading logarithmic terms. For example, the logarithm could be replaced 
by l/(6otts(9^))) and then the two-loop form of the running coupling could be used. 
These kind of changes should only produce small variations of the answer, if the scale 
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q is large enough for the leading logarithmic approximation to be working. If a large 
variation is found, this should be taken as an indication that subleading logarithmic 
corrections are large, and the result of the calculation should then only be taken as 
an estimate of the effect. 

The hadronic width of the r lepton has been used extensively for a determination 
of the strong coupling constant 113 . Since the r mass is only marginally large compared 
to typical hadronic scales, one should make sure that non-perturbative effects in r 
hadronic decays are under control. The r hadronic width is easily expressed in terms 
of the current- current correlators of axial and vector currents 



Rt- = Grri f dz{l — zY 

J\z\=l 



2ZMUP')+U^^lyip'] 



where z = p^/M^, and 



yrflU 



nSy(/)b>^ - p'gn + nSv'(/)p>'^ 



(6.1) 



(6.2) 



For the purpose of illustration we will now fix our reference values for the quark 
masses and condensates to those given in eq. ( ^.8|) . Defining 



Rr = i?l°) {l + Sj + <5f > + 5f > + 6\' 



m'' 



+ ... 



(6.3) 



(we do not include other QCD corrections in this formula) from formula ( |6.1|) and 



eq. (|5.7|) we obtain the following parametrization of the instanton corrections to Rr 



Si 



(boasirrir))-' 

-1 

11 
{boasirrir)) ^ 



rrir 



TTl-r 



im) 



5\' 



.{gg)' 



-1 



(6o«s("^r)) 

-1 

25" 

{boasirrir)) ^ 



l-39lg:o°lA3 
nir 



1.59tngA3 
rrir 



(6.4) 



where the numbers in parenthesis are only slightly sensitive to A3. Their central value 
was computed for A3 = 400 MeV, while the upper (lower) variations correspond to 
A3 = 500 MeV (A3 = 300 MeV). The range chosen for A3 corresponds to as{Mz) 
values of 0.113, 0.119 and 0.125, which is a reasonable representation of the present 
uncertainties. These results were obtained by replacing log(g^/A^) with l/(6oas(?^))) 
and then using the two-loop formula for the strong coupling constant in the complex 
plane integration. Numerical values for the corrections are given in table I, where the 
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Ag (MeV) 


300 


400 


500 


Si 


0.105 X lO-'^ 
(0.87 X 10-^) 


0.141 X 10-5 
(0.11 X 10-5) 


0.84 X 10-5 
(0.106 X 10-4) 


^m 


-0.99 X 10--^ 
(-0.60 X 10-5) 


-0.121 X 10-3 
(-0.70 X 10-^) 


-0.839 X 10-3 
(-0.46 X 10-3) 


^m^ 


-0.247 X 10-3 
(-0.10 X 10-3) 


-0.263 X 10-2 
(-0.95 X 10-3) 


-0.161 X 10-1 
(-0.51 X 10-2) 


^m^ 


-0.615 X 10-3 
(-0.14 X 10-3) 


-0.524 X 10-2 
(-0.82 X 10-3) 


-0.249 X 10-1 
(-0.16 X 10-2) 



Table 1: Instanton corrections to the r hadronic width 

values of the corrections obtained by using the one-loop form of the strong coupling 
constant are also shown in parenthesis for comparison. 

The corrections coming from the four- and six-quark condensates are by far the 
largest. They can reach the 2% level for the high end of the range of A3. This result 
is in qualitative agreement with ref. (we cannot expect exact agreement since 
we differ by subleading effects in our estimates) up to a sign, for which we have no 
explanation, since our analytical results do agree in sign. The four-quark condensate 
correction is of the same order as the six-quark correction. From these results we may 
conclude that instanton effects are not very important for the hadronic r decay, and 
that they are smaller than other sources of uncertainty. 

In ref. it was shown that instanton effects in the divergence of the correlator 
of two axial currents are large enough to spoil the standard methods to determine 
light quark masses from QCD some rulesB. The newly computed correction to the 
coefficient of the chiral condensate also contributes to these sum rules. The relevant 
quantity is 

^5(p')=PMP-n'^'- (6-5) 

With a straightforward generalization of the formalism of ref. H , defining 



m 



rriu + md 



rriu + rnd (. p 

^2^r^A 



33-2nf 



(6.6) 
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we obtain: 



TT 



,2_2 



27r 



11 
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p\9 rhs 1. 
a) q tt" 



52 /6.371AV/i3 




p 



p6 



log 



17 
P \ 9 



A 



(6.7) 



which extends the result of ref. 0. The contributions of the various terms of eq. ( ]67? 
to the finite energy sum rule is given by 

— cfc r/.s ImUfKf.sl = — 

27r 



— / ds Im^5(s) = A^o^' il + ^(*o) + ^(^o)} 

TT 1 J\s\=sn 47r 



s|=so 47r 

where R{so) is given in ref. M, and 



(6.8) 



^<^»)^fm':§hf' 



11 

9 



(6.9) 



It is easy to see that with fi = 180 MeV, T(so) is of order 1 for ^/sq ~ 4.1 A. For 
A = 400 MeV the instanton correction equals 1 already for sq ^ 2.7 GeV^ (while 
R{so) becomes 1 for sq ~ 2.2 GeV^). 

7. Conclusions 



In this paper we have presented a calculation of instanton corrections to the co- 
efficient function of chiral condensate operators in the operator product expansion 
of the correlator of fiavour non-diagonal axial or vector currents. The correction to 
the coefficient of the two-quark chiral condensate, together with the instanton correc- 
tion due to four-quark condensates in the factorization approximation is a new result 
of this work. We also repeated the calculation of the instanton correction due to 
six-quark condensates in the factorization approximation, and found agreement with 
ref. 1^. Our formula ( ^.7|) contains all the leading one instanton corrections to the 
correlator of off-diagonal vector or axial currents. 

The pattern of the corrections we found confirms the chiral nature of QCD. Thus, 
corrections proportional to the quark masses are smaller than corrections proportional 
to condensates, with the possible exception of corrections proportional to the strange 
quark mass and to the square of the chiral condensate, which may become larger than 
corrections proportional to the third power of the chiral condensate. 
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We applied our result to the calculation of the r hadronic width, where we found 
negligible corrections, and to the finite energy sum rules used to determine the light 
quark masses, where we found instead large corrections. Acknowledgements 
We would like to thank L. Girardello and G. Veneziano for useful conversations. 
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